Self-Consistent Vertex Correction Analysis for Iron-Based Superconductors: 
Mechanism of Coulomb-Interaction-Driven Orbital Fluctuations 
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We study the mechanism of orbital/spin fluctuations due to muhiorbital Coulomb interaction in 
iron-based superconductors, going beyond the random-phase-approximation. For this purpose, we 
develop a self-consistent vertex correction (SC-VC) method, and find that multiple orbital fluctua- 
tions in addition to spin fluctuations are mutually emphasized by the "multimode interference effect" 
described by the VC. Then, both the antiferro-orbital and ferro-orbital (=nematic) fluctuations si- 
multaneously develop for J/U ^ 0.1, both of which contribute to the s-wave superconductivity. 
Especially, the ferro-orbital fluctuations give the orthorhombic structure transition as weU as the 
softening of shear modulus Cee ■ 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20.Rp 
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Since the discovery of iron-based superconductors, the 
mechanism of high-Tc superconductivity has been studied 
very actively. Theoretically, both the spin-fluctuation- 
mediated s-t-wave state (with sign reversal of the gap 
between hole-pocket (h-pocket) and electron-pocket (e- 
pocket)) [1-5] and the orbital-fluctuation-mediated s++- 
wave state (without sign reversal) [6, 7] had been pro- 
posed. The latter scenario is supported by the robust- 
ness of Tc against impurities in many iron-pnictides [8- 
12]. Possibility of impurity-induced crossover from s± to 
S++ states had been discussed theoretically [3, 6]. Also, 
orbital-independent gap observed in BaFe2(As,P)2 and 
(K,Ba)Fe2As2 by laser ARPES measurement [13, 14] as 
well as the "resonance-like" hump structure in the neu- 
tron inelastic scattering [15] are consistent with the or- 
bital fluctuation scenario. 

Nature of orbital fluctuations has been studied inten- 
sively after the discovery of large softening of the shear 
modulus Cge [16-18] and renormalization of phonon ve- 
locity [19] observed well above the orthorhombic struc- 
ture transition temperature Ts- Consistently, a siz- 
able orbital polarization is observed in the orthorhombic 
phase [20, 21]. Moreover, the "electronic nematic state" 
with large in-plane anisotropy of resistivity or magneti- 
zation well above Ts and Tc [22, 23], also indicates the 
occurrence of (impurity- induced local) orbital order [24] . 

Origin of orbital order/fluctuation had been actively 
discussed, mainly based on the multiorbital Hubbard 
model with intra (inter) orbital interaction U (U') and 
the exchange interaction J = {U-U')/2 > [6, 25]. We 
had focused attention to a good inter-orbital nesting of 
the Fermi surfaces shown in Fig. 1 (a): Although moder- 
ate orbital fluctuations are induced by U' in the random- 
phase-approximation (RPA), the spin susceptibility due 
to the intra-orbital nesting, x'*(9)j is the most divergent 
for J > {i.e., U > U'). Since J/U « 0.12-0.15 accord- 
ing to the first- principle study [26], the RPA fails to ex- 
plain experimental "nonmagnetic" structure transition. 



This situation is unchanged even if the self-energy cor- 
rection is considered in the fluctuation-exchange (FLEX) 
approximation [27]. 

To explain the strong development of orbital fluctua- 
tions, we had introduced a quadrupole interaction [6]: 
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where g is the coupling constant, and O^ is the charge 
quadrupole operator ;7 = xz,yz,xy,x'^ — y^,3z'^ — r^. 
(Hereafter, a;,y-axes {X,Y-a.xes) are along the nearest 
Fe-Fe (Fe-As) direction.) This term is actually caused 
by the electron-phonon (e-ph) coupling due to in-plane 
Fe-ion oscillations [6, 14, 27]. Since O^ziyz) induces the 
inter-orbital scattering, strong antiferro (AF) orbital fluc- 
tuations develop for g > 0.2eV owing to a good inter- 
orbital nesting. We also studied the vertex correction 
(VC) beyond the RPA [28], and obtained strong en- 



hancement of ferro-quadrupole (O3. 



'J/Z 



susceptibility x^2_y2(0), which causes the orthorhom- 
bic structure transition and the softening of Cge [28]. 
This "nematic fluctuation" is derived from the interfer- 
ence of two AF orbitons due to the symmetry relation 
d^2_y2{0) ^ OxziQ) X Oyzi-Q), where OxziYZ) = 
[Oxz + {~)Oyz]/V^- Then, it was natural to expect 
that such multi-orbiton interference effect, which is given 
by the VC while dropped in the RPA, induces large 
"Coulomb-interaction-driven" orbital fluctuations. 

In this letter, we study the orbital and spin fluctua- 
tions in iron-based superconductors by considering the 
multiorbital Coulomb interaction with U ^ U' + 2J and 
J/U - 0(0. 1). We develop the self-consistent- VC (SC- 
VC) method, and find that both ferro-Oj,2_y2 and AF- 
Oxz/yz fluctuations strongly develop even for J/U ^0.1, 
due to the inter-orbital nesting and the positive in- 
terference between multi-fluctuation (orbiton-f magnon) 
modes. This result leads to a conclusion that RPA un- 
derestimates the orbital fluctuations in multiorbital sys- 



terns. The present study offers a unified explanation for 
both the superconductivity and structure transition in 
many compounds. 

Here, we study the five-orbital Hubbard model intro- 
duced in Ref. [1]. We denote d-orbitals m = 3z^ — r^, 
xz, yz, xy, and x^ — y^ as 1, 2, 3, 4 and 5, respectively. 
The Fermi surfaces are mainly composed of orbitals 2, 3 
and 4 [28]. Then, the susceptibility for the charge (spin) 
channel is given by the following 25 x 25 matrix form in 
the orbital basis: 



r^'\q) = r'^''^'\q){i - t<''^r''^''^'\q)r\ 



(2) 



where q = {q,^i = 2ttIT), and F'^^'') represents the 
Coulomb interaction for the charge (spin) channel com- 
posed of C/, U' and J given in Refs. [6, 14]. The irre- 
ducible susceptibility in Eq. (2) is given as 



r''<^\q) = x\q)+X<^\q). 



(3) 



where Xu^mm'iq) = ~TY.pGi^{p + q)G„,'i'{p) is the 
bare bubble, and the second term is the VC (or orbiton 
or magnon self-energy) that is neglected in both RPA 
and FLEX approximation. In the present discussion, it 
is convenient to consider the quadrupole susceptibilities: 



IV ,mm' 

= Tr{0^r((?)Oy}. 
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Non-zero matrix elements of the quadrupole operators 



for the orbital 2 
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1 [28]. Because of the symmetry, the off- 
diagonal susceptibilities (7 ^ 7') are zero or very small 
for q = and the nesting vector Q « (tt, 0) or Q' w (0, it) 
[28]. We do not discuss the angular momentum (dipole) 
susceptibility, x'jj.iQ) "^ (^m(9)^m ("?))' since it is found to 
be suppressed by the VC. Note that O^^ oc l^l,y + 1^1 ^i- 

To measure the distance from the criticality, we in- 
troduce the charge (spin) Stoner factor Uq , which 
is the largest eigenvalue of r'^'^'*)x"''''°(*')(q') at uji — 0: 
The charge (spin) susceptibility diverges when ctmax = 
maxq{aq } = 1. In a special case J = 0, the relation 
Q^max = c«max holds at the momentum Q in the RPA; see 
Fig. 1 (b). That is, both spin and orbital susceptibilities 
are equally enhanced at J = 0, which is unchanged by the 
self-energy correction in the FLEX approximation [27]. 
For J > 0, the spin fluctuations are always dominant 
(•^max > Q^max) i^i the RPA Or FLEX. However, because 
of large X'^{q), the opposite relation a^^x ^ o^max ^^"^ be 
realized even for J/U < 0.1 in the SC-VC method. 

First, we perform the RPA calculation for n — 6.1 and 
T — 0.05, using 32 x 32 fc-meshes: The unit of energy is 
eV hereafter. Figure 1 (c) shows the diagonal quadrupole 
susceptibilities for J/U = 0.088; x^iq) = X'^-yiq)- (The 
spin susceptibility is shown in Ref. [1].) The Stoner 
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FIG. 1: (color online) (a) Fermi surfaces of iron pnictides. 
The colors correspond to 2 — xz (green), 3 — yz (red), and 
4 = xy (blue), respectively, (b) qq, qq and U as function of 
J/U in RPA under the condition ofnax = 0.97. (c) Xxzil) and 
Xl^-y2{q) given by the RPA for {J/U, U) = (0.088, 1.53). 



factors are a^ax = 0.97, Uq = 0.76, and Uq = 0.47; 
see Fig. 1 (b). In the RPA, xlziQ) IxlziQ')] is weakly 
enlarged by the intcr-orbital (3,4) [(2,4)] nesting, while 
Xx2-i;2 (9) i^ relatively small and AF-like. Thus, the RPA 
cannot explain the structure transition that requires the 
divergence of x'x2^y2{0)- 

Next, we study the role of VC due to the Maki- 
Thompson (MT) and Aslamazov-Larkin (AL) terms in 
Fig. 2 (a), which become important near the critical 
point [29, 30]. Here, X<'\q) = X^^'^q) + (-)X'^'-^(q), 
and wavy lines represent x"'^- The AL term (AL1-I-AL2) 



for the charge sector, X 
T 



AL,c 
W .mm 



(q), is given as 



+3VXcdik + q)V:f^ghi-k)}K^„,,^,a^g,,{q; k), (5) 

where V^^^q) = f^.= + f'*^=x'*'''('?)^"'^ A((j; fc) is the 
three-point vertex made of three Green functions in Fig. 
2 (a) [28], and Kim',cd,ghiq'^k) = Kch,mg,dm'{q;k) + 
^gd,mcjim'{q', —k — q). We includc all [/^-terms, which are 
important for reliable results. The expressions of other 
VCs will be published in future. 

Both MT and AL terms correspond to the first-order 
mode-coupling corrections to the RPA susceptibility: 
The intra- (inter-) bubble correction gives the MT (AL) 
term [29]. In single-orbital models, the VC due to 
MT-I-AL terms had been studied by the self-consistent- 
renormalization (SCR) theory [29] or FLEX approxima- 
tion with VC [30], and successful results had been ob- 
tained. In the former (latter) theory, the susceptibility 
is calculated in the self-consistent (self-inconsistent) way. 
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FIG. 2: (color online) (a) The MT and AL terms: The wavy 
and solid lines are susceptibilities and electron Green func- 
tions, respectively. A;;',a6,e/ is the three-point vertex, (b) 
Dominant AL terms for Xi2_„2(0); the first (second) term 
represents the two-orbiton (two-magnon) process, (c) Domi- 
nant AL terms for xliiQ) (-^ = 2:2, j/z); higher-order terms 
with bubbles made of Xm (±Q) (= multi-fluctuation process) 
are relevant, (d) Enhancement of An',a6,e/ due to charge VCs. 



Here, we find a significant role of the AL term inherent 
in the multiorbital Hubbard model. 

Now, we perform the SC-VC analysis, in the way to 
satisfy x'^'^'iq) in the VC are equal to the total suscepti- 
bilities in Eq. (2). Then, x'^{q) is strongly enhanced by 
j^AL.c jj^ gq^ ^g^^ which is relevant when either x'^ or x'* 
is large. On the other hand, we have verified numerically 
that X" - TX; A • VW ■ A is less important, although 
it could be relevant only when both x'^ and x'' are large. 
Hereafter, we drop X'^{q) to simplify the argument. Fig- 
ure 3 (a) show x'^{q) given by the SC-VC method for 
n = 6.1, J/U = 0.088 and U = 1.53, in which the Stoner 
factors are a^ax = 0^0 = 0.97 and ag = 0.86. Compared 
to the RPA, both Xx-^-y^io) and x'xzi'i) are strongly en- 
hanced by the charge AL term, X^^^'^, since the results 
are essentially unchanged even if MT term is dropped. 
In the SC-VC method, the enhancements of other charge 
multipole susceptibilities are small. Especially, both the 
density and dipole susceptibilities, ^i mXii mmi'i) ^^^ 
xliil) (m — 2;, 2/, 2) respectively, are suppressed. 

Here, we discuss the importance of the AL 
term: At q ~ or Q, x'^io) is enlarged by 
the diagonal vertex correction with respect to 7, 
X^^'^iq) = Tr{6^XAL,c(q)0^}/Tr{62}, since the off- 
diagonal terms are absent or small [28]. The charge AL 
term in Eq. (5) is given by the products of two x'^'s 
(two-orbiton process) and two x'''^ (two-magnon pro- 
cess), shown in Fig. 2 (b). The former process was 



discussed in Ref. [28], and the latter has a similar- 
ity to the spin nematic theory in Ref. [16] based on 
a frustrated spin model. Now, we consider the orbital 
selection rule for the two-orbiton process: Because of 
the relation Tr{Oa;2_j,2 0^j} ^ for M — xz,yz and 
a rough relation Kw ,ab,cd ^ ^w MbM^v ,a5b,c6d,i [28], the 
two-orbiton process for ^ — x^ — y^ is mainly given 
by xliiQf- According to Eq. (5) and Ref. [28], 
^^2^42(0) - A2[/4T^^{x(g)}2 grows in proportion to 
TxiQ) [log{x(Q)}^] at high [low] temperatures. In the 
case of Fig. 3 (a), two-magnon process is more impor- 
tant for x^2_„2(0) because of the relation Uq > Uq. We 
checked that the two-magnon process is mainly caused 

byxi2,22(Q)'-xi2,33(Q)'>o. 
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FIG. 3: (color online) XxzIq) and Xx^-y^iQ) given by the SC- 
VC method. The relation a^ax = ceo = 0.97 is satisfied in all 
cases: (a) n = 6.1 and J/U = 0.88 (qq = 0.86), (b) n = 6.1, 
J/U = 0.9 and g = 0.05 (oq = 0.96), and (c) n = 6.05, 
J/U = 0.11 and g = 0.065 (a^ = 0.87). 

In the same way, X£,(Q) - ^^U^TJ2^XMil + 
Q)x'x2_y2{q) is enlarged by the two-orbiton process due 
to xliiQ) and x^2_j^2(0), shown in Fig. 2 (c). (In this 
case, two-magnon process is less important since x^(0) 
is small.) The obtained Xxz{q) has peaks a.t q = Q 
and Q' since the inter-orbital scattering is emphasized 
by X^,{Q) (X xS2_,2(0) » 1. Thus, both XxAQ) and 
Xx2-y2{0) are strongly enlarged in the SC-VC method, 
because of the "positive feedback" brought by these two 



AL terms: Figure 2 (c) shows an example of the higher- 
order terms that are automatically generated in the SC- 
VC method. Such "multi-fluctuation processes" inherent 
in the self-consistent method magnify the RPA results. 

Thus, strong ferro- and AF-orbital fluctuations are 
caused by AL terms. Both fluctuations work as the 
pairing interaction for the s_|__|_-state, while the fcrro- 
fluctuations are also favorable for the s±-state. For 
J/U < iJ/U)c = 0.088, the relation a^^^ < ^o = 0.97 
is realized and Oq increases towards unity. In this case, 
orbital order occurs prior to the spin order as increasing 
U with J/U is fixed, since the VC (due to two-orbiton 
process) can efficiently enlarge orbital susceptibilities be- 
cause of large a^ax(I^PA)- This situation would be con- 
sistent with wider non-magnetic orthorhombic phase in 
Nd(Fe,Co)As and many 1111 compounds. 

Since the present SC-VC method is very time- 
consuming, we applied some simplifications: We 
have verified in the self-inconsistent calculation that 
Ti{0~^X{q)0~^'} with 7 7^ 7' is zero or very small, 
especially at q = and Q for the reason of sym- 
we are interested in the enhance- 
at qr = and Q and the dom- 



between 7 = xz,yz,x'^ — y^ , 



nictry. Since 

ment of x^(q) 

inant interferences 

we calculated ^;/',mm'('?) only for {(Z, Z'), (m,m')} e 

xz,yz,x^ — y"^. [{ij') G 7 means that Ol^'^ ^ 

0.] That is, {(;,n,Km')} = {(1, 2), (3,4), (2, 5)}, 

{(1, 3), (2, 4), (3, 5)}, and {(1, 5), (2, 2), (3, 3)}. 

We stress that both {J/U)c and AF-orbital fluctua- 
tions increase by considering following two factors: The 
first one is the charge VC at each point of the three- 
point vertex in Fig. 2 (d), as a consequence of the Ward 
identity between A and x^". The enhancement factor 
at each point is estimated as 1 + XH^/xj = 1.3 ^ 2.5 
for J = xz and x'^ — y^ in the present calculation near 
the critical point. This effect will increase {J/U)c sen- 
sitively. The second factor is the e-ph interaction: We 
introduce the quadrupole interaction in Eq. (1) due 
to Fe-ion oscillations [6, 14, 27]. As shown in Fig. 3 
(b), very strong AF-orbital fluctuations are obtained for 
J/U = 0.09 and g = 0.05; a^^^ = ajj = 0.97 and 
Uq — 0.96. The corresponding dimensionless coupling 
is just A = gN{0) - 0.035 [6, 27]. We also study the 
case n = 6.05 and g — 0.065, and find that the relation 
Q^max = Q^max = ^■^'^ '^^ realized at {J/U)c = 0.11, as 
shown in Fig. 3 (c). For these reasons, strong ferro- and 
AF-orbital-fluctuations would be realized by the cooper- 
ation of the Coulomb and weak e-ph interactions. 

Finally, we make some comments: The present multi- 
fluctuation mechanism is not described by the dynamical- 
mean- field theory (DMFT), since the irreducible VC is 
treated as local. Also, the local density approximation 
(LDA), in which the VC is neglected, does not reproduce 
the nonmagnetic orthorhombic phase. Although Yanagi 
et al. studied U' > U model [7] based on the RPA, 
that was first studied in Ref. [31], x%z^-r^i^) develops 



while Xs2_w2(0) remains small, inconsistently with the 
structure transition. Our important future issue is to 
include the electron self-energy correction into the SC- 
VC method, which is important to discuss the filling and 
T-dependences of orbital and spin fluctuations, and to 
obtain more reliable {J/U)c- 

In summary, we developed the SC-VC method, and 
obtained the Coulomb-interaction-driven ncmatic and 
AF-orbital fluctuations due to the multimodc (or- 
bitons-|-magnons) interference effect [28] that is over- 
looked in the RPA. For J/U < {J/U)c, the structure 
transition [a^ « 1) occurs prior to the magnetic tran- 
sition [uq Ri 1), consistently with experiments. When 
Q^max ^ '^max: both S++- and s±-states could be realized, 
depending on model parameters like the impurity con- 
centration [3,6]. In a sense of the renormalization group 
scheme, the quadrupole interaction in Eq. (1) is induced 
by the Coulomb interaction beyond the RPA. We expect 
that orbital-fluctuation-mediated superconductivity and 
structure transition are realized in many iron-based su- 
perconductors due to the cooperation of the Coulomb 
and e-ph interactions. 

This study has been supported by Grants-in-Aid for 
Scientific Research from MEXT of Japan, and by JST, 
TRIP. Part of numerical calculations were performed on 
the Yukawa Institute Computer Facility. 
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